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Abstract 

We perform a comparative study of two methods of determining the survival prob- 
abilities of low, intermediate, and high energy solar neutrinos that emphasizes the 
general agreement between the Large Mixing Angle (LMA) solution and extant solar 
neutrino data. The first analysis is oscillation parameter-independent and the second 
analysis involves an approximate calculation of the survival probabilities in the three 
energy ranges that depends only on oscillation parameters. We show that future ex- 
periments like BOREXino, CLEAN, Heron, LENS and MOON, that measure pp and 
7 Be neutrinos, will facilitate a stringent test of the LMA solution independently of the 
Standard Solar Model (SSM), without recourse to earth-matter effects. Throughout, 
we describe the role of SSM assumptions on our results. If the LMA solution passes the 
test without needing to be modified, it may be possible to establish that 9 X is nonzero 
at more than 2a assuming the SSM prediction for the pp flux is correct. 
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1 Introduction 



The LMA solution to the solar neutrino problem has emerged as the preferred explanation 
of solar neutrino data PUI2|- Reactor neutrino data from the KamLAND experiment |3, have 
lent much confidence in this solution Nevertheless, it has long been recognized that the 
LMA solution is not in wholly satisfactory agreement with solar neutrino data pp. Since 
the media traversed by neutrinos incident at KamLAND and at solar neutrino experiments 
are considerably different, it is important to confirm that the LMA solution with matter 
effects dictated by the MSW mechanism is consistent with solar neutrino data. Recent work 
emphasizing our ignorance of neutrino-matter interactions and suggestions on how solar data 
can help illuminate the nature of these interactions can be found in Ref. [3]. 

For the LMA solution, the Chlorine measurement with its non- negligible component 
of 7 Be neutrinos, is too low to be consistent with the SuperKamiokande [7j (SuperK) and 
SNO [H| measurements of the 8 B neutrinos. The Small Mixing Angle solution can accommo- 
date a low survival probability of 7 Be neutrinos, but at the expense of a highly non-uniform 
suppression of the 8 B neutrinos relative to the SSM [D] spectrum, which is contradicted by 
SuperK and SNO data. Exotic mechanisms have been proposed to account for the suppres- 
sion of the 7 Be neutrinos, but these require that either a sterile neutrino be added [TU] to 
the standard three- neutrino framework or that neutrino masses vary with density [TTj . 

In this Letter we perform a model-independent analysis of the latest solar neutrino data 
with the flux normalizations and survival probabilities of the low, intermediate, and high en- 
ergy neutrinos all treated as independent parameters; the possibility of a model-independent 
test of the SSM was emphasized early on in Ref. ^2]- We then determine which of these pa- 
rameters are calculable from current data and which require additional input from the SSM. 
The results of this analysis are then used to test the LMA predictions for neutrinos in each 
energy range. We show that with the predictions of the latest SSM and including the latest 
SNO salt data ^3] the agreement between the LMA solution and the solar data remains 
relatively good, although there are minor discrepancies between the oscillation probabilities 
of the low- and intermediate-energy solar neutrinos which cannot be resolved by extending 
the analysis to include the mixing of a third neutrino. We demonstrate that future solar 
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neutrino experiments that would measure low and intermediate energy neutrinos and reduce 
the uncertainties in the survival probabilities will provide a critical test of the LMA solution. 



2 Current data 

2.1 Model independent analysis 

2.1.1 Formalism 

We use the following notation: R is the ratio of the measured rate to the SSM prediction 
for a given experiment, f3 is a flux normalization relative to the SSM, and Pl, Pi and Ph 
are average survival probabilities of low energy (pp), intermediate energy ( 7 Be, pep, 15 0, 
13 N) and high energy ( 8 B, hep) neutrinos, respectively. For three active neutrinos and the 
recently updated SSM 0, the relative rates are given by 

R Ga = 0.109/3^ + 0.335/3^ + 0.556/3^, (1) 

Rci = 0.803/3^ + 0.197/3^, (2) 

Rsno = PhPh , (3) 

R N sSo = Ph- (4) 

Here, CC and NC refer to charged- current and neutral-current measurements, respectively. 
See Refs. |14[ ITS] for a description of the method. The present measurements of these 
quantities are given in Table HJ The elastic scattering rates at SuperK and SNO are not 
included since they provide redundant information with less precision than Rg^o anc ^ Rsnoi 
although in principle these data could also be included to improve the accuracy. 

2.1.2 Without SSM constraints 

The quantities PlPl, PiPi and PhPh are determined from the CC measurements RGa, Rci 
and Rsno- Rsno determines Ph- Using the data in Table HI we find 

(3 L P L = 0.69 + 0.11, (5) 
/3/P/ = 0.37 + 0.15, (6) 
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Table 1: Current solar neutrino measurements. 
Measurement value source 

R Ga 0.54 ± 0.03 SAGE/GALLEX/GNO [IE] 

Rci 0.31 ± 0.03 Homestake |6! 

Rgfio 0.30 ± 0.01 SNO D 2 phase, SNO salt phase HSJ 
r sno °- 88 ± °- 06 SN0 D 20 phase, SNO salt phase EE2| 

p H p H = 0.30 ±0.01, (7) 
fa = 0.88 ±0.06, (8) 



where the uncertainties are la. Then 



Ph = RTno/Rsno = 0-34 ± 0.03 . (9) 

We note that Ph would not be determined by Rsno an d Rsno a l° ne if there were a fourth, 
sterile neutrino since the expression for R^no w °uld depend on the sterile content of the 
oscillating neutrinos [T^]. At present there are no data that can determine fti and Pi 
(or (3i and Pi) separately in a completely model-independent way without imposing SSM 
constraints. 



2.1.3 With SSM constraints 

The la fractional uncertainties of the solar neutrino fluxes from the SSM jHj are 

5{P S L SM ) = 0.010, (10) 

5{j3f SM ) = 0.122, (11) 

5{p H SM ) = 0.163. (12) 

The values of Pl and Pj can be determined from current data only if the SSM constraints 
are used; from Eqs. (0 El HSJ) an d (fTTj) we deduce 

P L = 0.69 ±0.11 with SSM flux, (13) 

Pi = 0.37 ±0.16 with SSM flux. (14) 
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If the SSM prediction for (3jj is used as an additional constraint, then a combination of 
the SSM and Rgg gives 

Ph = 0-89 ± 0.06 , (15) 



from which we deduce 



P H = 0.34 ± 0.03 . (16) 



2.2 LMA analysis 
2.2.1 Formalism 

If the v e disappearance of solar neutrinos occurs via the MSW mechanism [T7| with adiabatic 
propagation jTHJ (which is the situation for the LMA solution), the solar neutrino oscillation 
probability in a three-neutrino framework is given by the approximate formula [TH] 1 . 



ci 



1 + cos 20™ cos20 s 
2 



+ at . (17) 



The above equation applies when 9 X is small, as is indicated by the CHOOZ experiment |21j ; 
the 95% C. L. is sin 2 ^ < 0.05 for the best-fit atmospheric mass-squared difference of 0.0021 
eV 2 from the SuperK experiment [22 a . For the current range of LMA parameters preferred 
by data, earth-matter effects [23] are very small |Hj and therefore we do not consider them 
here. The quantity in the square brackets is the two-neutrino oscillation probability and the 
factors involving 9 X are the corrections due to mixing with the third neutrino. The angle 9™ 
is the mixing angle in matter at the point of origin of the neutrino, and is given by 

tan 20?= 8in29s (18) 
cos29,-A 



2V2G F N°E u cl _ t n ^ . / E v \/8x lO" 5 eV 2 \ / N° e . , 



where 

As " ~~Lr~* - 19 * 10-3 (t^v ){ .v„/ j ) * • < ig > 

The initial electron densities at the neutrino source are approximately N® /(iV^/cm 3 ) = 106, 

86 and 58 for the high, intermediate and low energy neutrinos, for which the corresponding 
^dn our notation, 5m 2 s is the solar mass-squared difference and 9 S and 9 X are the mixing angles conven- 
tionally denoted by #12 and #13, respectively |2()j . 
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average neutrino energies are 9.06, 0.862, and 0.325 MeV, respectively. The factor c 2 in 
A is a three-neutrino correction to the effective electron number density [T""]. In the LMA 
analysis with three neutrinos, the variables that determine the probabilities are therefore 
5m 2 , 9 S and 9 X , instead of the oscillation parameter-independent probabilities Pj. 

2.2.2 Without SSM constraints 

As discussed in the model-independent analysis, neutrino data alone cannot presently de- 
termine Pl or Pj, so these probabilities cannot be used to constrain the LMA parameters. 
However, the SNO CC and NC data can be used to determine Ph, and 9 S can then be 
determined via Eqs. ()17H19j) if a range of 5m 2 is used as input and we ignore 8 X . Taking the 
la range, 5m 2 = (7.9 ± 0.55) x 10" 5 eV 2 , from KamLAND data [3J, we find 

tan 2 9 S = 0.45 ±0.06, (20) 

in good agreement with recent global fits to 9 S that also include day/night spectral shapes 
in the analysis [21]. If 9 X is allowed to be nonzero, 9 S can be calculated as a function of 8 X ; 
see Fig. An increasing 9 S can be compensated by an increasing 9 X . 

Since fti and f3i are as yet undetermined by the data, we can use the LMA prediction for 
Pl and Pi in conjunction with the current best-fit values of /3lPl and fiiPi from Eqs. (J""} 
and (0) to solve for (3l and (3i (see Fig. [2]); the LMA predictions assume the la ranges, 
Sm 2 s = (7.9 ± 0.55) x 10~ 5 eV 2 from KamLAND and tan 2 9 S = 0.45 ± 0.06 from Eq. (^U|) . 
We see that the flux normalization of the low (intermediate) energy neutrinos needs to be 
slightly higher (lower) in order for the LMA predictions to be completely consistent with 
current data, although the present uncertainties are sufficiently large that there is essentially 
no conflict. 

2.2.3 With SSM constraints 

The LMA predictions for the solar neutrino survival probability versus neutrino energy 
are shown in Fig. [3] for sin 2 #3. = and 0.05 with 5m 2 and 9 S varying over their currently 
preferred ranges. Current data (taken from Eqs. 19*1 fTBI and [T4*|) are also shown, where the SSM 
constraints in Eqs. ("!"""") and ("IT]) have been used to determine Pl and Pi. The oscillation 
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probability of the low-energy (intermediate-energy) neutrinos is slightly higher (lower) than 
the LMA prediction, although the discrepancies are currently only at the la level. Increasing 
9 X lowers the LMA curve, which improves the fit to the intermediate energy data, but at 
the expense of a worse fit to the low energy data. Therefore reductions in the uncertainties 
from future solar neutrino experiments will provide a critical test of the LMA solution (see 
Sec. 3). 

Using the LMA expressions for P^, Pj and Ph, contours of constant <5m 2 and tan 2 9 S are 
plotted in Fig. HI in the (Ph, Pl) and (Ph, Pi) planes for sin 2 9 X = and 0.05. The data 
points and larger error bars are taken from Eqs. (0 IT3*|) and (|14j) . and the smaller error 
bars illustrate the expected uncertainties from future measurements of Rp P c , Rp P S , R%ei 
-Rff, Rsno an d R-sno whh SSM constraints imposed (see the next section). The effect of 
increasing 9 X on the compatability of the LMA solution with the data is evident. 



3 Future data 

3.1 Model independent analysis 
3.1.1 Without SSM constraints 

Future experiments such as MOON [25 or LENS j2H] ot will be able to provide better 
measurements of CC scattering of the low and intermediate energy neutrinos, i.e., of 

Rpp — PlPl, (21) 

RBe = PlPl, (22) 

perhaps at the 2.5% level [2Zj. However, these measurements will still not separate the 
flux normalization from the survival probability. In order to determine the flux normaliza- 
tion from data alone (i.e., without imposing theoretical inputs from the SSM), a process 
with a NC component must be used, such as the elastic scattering (ES) measurements at 
BOREXino [28 , CLEAN [H], HERON [H] and KamLAND 31 . They will measure 

R^ = (3 L P L (l-r L )+(3 L r L , (23) 
Rfe = PiPiQ-rfi + Pirj, (24) 
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Table 2: Expected uncertainties from future solar neutrino measurements of CC, NC and 
ES processes. 





pp 


7 Be 8 B 


CC 


2.6% 


2.5% 3.5% 


NC 




- 4.8% 


ES 


1.1% 


2.6% - 



where rj is the ratio of the NC to v e cross-sections for scattering off electrons in the appropri- 
ate energy regime. CLEAN plans to measure i?£f to 1.1% and R^f to 2.6%, HERON plans 
to measure R pp s to 3.5%, and BOREXino plans to measure R B ^ to 5% |27j . KamLAND can 
measure i?ff to about 11% [SI]. Once R pp s and are measured, all six parameters (three 
flux normalizations and three survival probabilities) will be determined by neutrino data. 
Furthermore, SNO expects to reduce the uncertainties on its CC and NC measurements to 
about 5.5% and 6.4%, respectively, in its third phase J22]- When combined with its previous 
measurements, the SNO uncertainties will be about 3.5% and 4.8%, respectively. The future 
expectations for these uncertainties are summarized in Table El where we list the smallest 
uncertainty in each channel. 

The pp measurements can be used to determine Pl and Pl via 

p L = [Rg s -R™(i-n)]/n, (25) 

Pl = flgprj/ [R E J - - r f )] , (26) 

and similarly Pi and Pj can be determined from R B ^ and R B ^. Using the six measurements 
Rpp, Rsei P-snoi P-snoi Ppp an d P-Be with the projected uncertainties shown in Tabled the 
values of Pl, Pi and Ph can be determined independently from any solar model assumptions, 
with uncertainties of about 16%, 22% and 5%, respectively, assuming the central values of 
the PjPj and Ph remain about the same (Eqs. EHHJ) an d that the best-fit values of Pl and Pi 
are close to unity. The constraints on Pl and Pi are not nearly as precise as the pp and 7 Be 
measurements themselves because the NC component, from which the value of P is inferred, 
is suppressed by the smaller NC cross-section for and z/ T , compared to the ES cross-section 
for u e . The resulting uncertainties of Pl, Pi and Ph would be about 14%, 12% and 6%, 
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respectively, if solar model constraints are not imposed. 
3.1.2 With SSM constraints 

In the future, using the CC measurements of the pp and 7 Be neutrinos to better determine 
PlPl and /3/Pj, respectively, the uncertainties of Pi and Pj will be reduced to 3% and 11%, 
respectively, when the SSM constraints are used; the uncertainty of Pj does not improve much 
since the SSM intermediate energy flux is known to only 12%. The large SSM uncertainty 
of Ph in Eq. ()12j) does not provide any appreciable reduction in the overall uncertainty of 
(3h when combined with SNO NC data. 

3.2 LMA analysis 

3.2.1 Without SSM constraints 

Once the six measurements in Tableware made, the six parameters /3l, /3j, (3h, Sm^, 9 S and 
9 X will in principle be determined from Rs^ , Rsnoi ^Be ' ■> -^sfj Rpp an d Rpp > without 
needing any other inputs such as SSM fluxes or KamLAND data. 

Figure El shows contours of constant 9 S and 9 X in two-dimensional subspaces of probability 
space for 5m 2 s = 8 x 10~ 5 eV 2 . The current best fit values for these parameters are also 
shown in the figure. The central value of Pi is accommodated only for sin 2 ^. < 0, which is 
unphysical, while the central value of Pi is accommodated for values of sin 2 9 X that violate the 
CHOOZ bound. However, there is adequate consistency within error bars. As noted earlier, 
if the flux normalizations (3l and /3j of the SSM are correct, then the current central values 
of Pl and Pj are slightly high and low, respectively. Future solar neutrino measurements will 
be able to provide a more definitive statement about this minor discrepancy, even without 
making any SSM assumptions. It is also clear from Fig. |3] that a more precise experimental 
determination of (3j and especially f3i would be needed to provide significant improvements 
in the determination of oscillation parameters if the predictions of the SSM are not used. 
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3.2.2 With SSM constraints 

In the future, after is measured and the SSM constraint (Eq. [TUJ) is imposed, Pi would 
be determined to about 3%. These projected uncertainties are shown as the smaller error 
bars in Figs. 0] and 03 

With the SSM constraints imposed, measurements of tan 2 9 S and sin 2 9 X can be made with 
la uncertainties of order 0.03 and 0.015, respectively, by future solar neutrino experiments. 
As is evident from Figs. 0] and to be consistent with the LMA solution, the current central 
values of Pi and Pi should shift when the new measurements are made. 

We note that if Rci were about la higher, then both PlPl and (3iPj would be about at 
the values predicted by the best-fit LMA parameters and the SSM; they are coupled together 
since PlPl cannot be determined from Rc a without knowing Pi Pi, which is determined from 
Rci and PhPh = Rsno- Since the primary constraint on Pi Pi comes from R C i, to which fiiPi 
only makes a 20% contribution, the current determination of (3iPi has a large uncertainty. 
Future measurements of (3iPi in R^ will not only significantly reduce the uncertainty on 
Pi Pi, but also will test the Chlorine measurement. 

4 Summary 

We have shown that current neutrino data can determine the survival probability of only 
the high-energy neutrinos if SSM flux normalizations are not used, and that if the SSM 
flux constraints are imposed, the LMA solution provides a consistent explanation of all 
solar neutrino data, although minor discrepancies remain at the la level. In particular, the 
implied flux of low (intermediate) energy neutrinos appears to be slightly high (low), if the 
LMA solution is correct. Including the small mixing effects of the third neutrino cannot 
compensate for both inconsistencies simultaneously. 

Future measurements of low and intermediate energy neutrinos will provide a much 
stricter test of the viability of the LMA solution independently of SSM predictions for the 
fluxes. See Table |3]for a summary of our results. If reliance is placed in the SSM pp flux nor- 
malization, the test becomes conclusive. If the LMA solution survives without modification, 
it may be possible to demonstrate that 9 X is nonzero at more than 2a. 
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Data 


SSM imposed? 


<W 




<w 


<W 




S(Ph) 


Current 


No 






6 






8 


Current 


Yes 


1 


12 


6 


15 


43 


8 


Future 


No 


16 


22 


5 


14 


12 


6 


Future 


Yes 


1 


11 


5 


3 


11 


6 



Table 3: Percentage uncertainties in the flux normalizations and survival probabilities from 
current and future solar neutrino measurements. From the current data, Pl and Pi (and 
hence Pl and Pj) cannot be determined without SSM input. For the future data we have 
assumed that the central values of /3jPj and Ph coincide with those of the current measure- 
ments (see Eqs. [3HHJ) and that the best-fit values of Pl an d Pi are close to unity. 
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Figure 1: Inferred values of tan 2 9 S versus sin 2 9 X using the current SNO data and the la 
range, 5m 2 s = (7.9 ± 0.55) x 10" 5 eV 2 , from KamLAND data alone. 
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Figure 2: Inferred values of (3l (solid) and /?/ (dashed) for 5m 2 s = (7.9 ± 0.55) x 10~ 5 eV 2 , 
using the LMA predictions for Pl and Pi and Eqs. (|5I6|) . In each case the central values and 
la error band is shown. The dotted line is the SSM prediction. 
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Figure 3: LMA predictions for solar neutrino survival probability versus neutrino energy 
for sin 2 6 X = (solid), 0.05 (dashed). The bands shows the range of predictions for 8m 2 s = 
(7.9 ± 0.55) x 10~ 5 eV 2 and tan 2 6 S = 0.45 ± 0.06. The data points are taken from Eqs. (0 
HHJ) and (JHJ), the latter two of which were obtained using the SSM predictions for f3 L and 
Pi. 
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Figure 4: Contours of 9 S and 5m 2 s in the (Ph,Pl) an d (Ph,Pi) planes for sin 2 ^ = 
and 0.05. The bold data point and larger error bars are taken from Eqs. (0 fT^j) and (|14|) . 
the latter two of which were obtained using SSM constraints. The smaller error bars are 
projected improvements in the uncertainties from future measurements of Rpp, R%ei 
Rge, Rsno an d Rsno w ^h SSM constraints imposed. 
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Figure 5: Contours of 9 S and 9 X in the (Ph, Pl) and (Ph, Pi) planes, for 5m 2 s = 8 x 10~ 5 eV 2 . 
The bold data points are the current central values from Eqs. (JHJ IT3*j) and ()14j). The larger 
(smaller) error bars are the projected uncertainties from the future measurements listed 
in Table 121 without (with) SSM constraints imposed. While adding SSM constraints to 
future data reduces the uncertainty in Pl significantly, Pi and Ph are essentially unim- 
proved. Two representative future measurements (with SSM flux constraints) are indicated 
at (tan 2 9 S , sin 2 X ) = (0.35,0.03) and (0.44,0). 
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